Transport of free surface liquid films and drops by external ratchets and self- ratcheting mechanisms 
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We discuss the usage of ratchet mechanisms to transport a continuous phase in several micro-fluidic settings. 
In particular, we study the transport of a dielectric liquid in a heterogeneous ratchet capacitor that is periodically 
switched on and off. The second system consists of drops on a solid substrate that are transported by different 
types of harmonic substrate vibrations. We argue that the latter can be seen as a self-ratcheting process and 
discuss analogies between the employed class of thin film equations and Fokker-Planck equations for transport 
of discrete objects in a 'particle ratchet'. 
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I. INTRODUCTION 

Most of the many studied ratchet mechanisms are discussed 
in the context of directed transport and filtering of discrete 
objects as, e.g., colloids or macromolecules llll. Examples 
include colloidal particles, that move in a directed manner 
through an array of periodic asymmetric micropores when un- 
der the influence of an oscillating external pressure Such 
particles may as well be driven by a dielectric potential of 
sawtooth shape that is switched on and off periodically ^ . In 
many of those 'discrete' ratchets the carrier fluid (or solvent) 
does not or nearly not move on average, i.e., they can not be 
employed to transport a continuous phase. One exception is 
the motion of magnetic particles in ferrofluids under the influ- 
ence of an oscillating magnetic field There the resulting 
net motion of the particles is transmitted to the carrier liquid 
by a strong viscous coupling. 

However, the basic concept of ratchet transport - that a lo- 
cally asymmetric but globally homogeneous system may in- 
duce global transport if it is kept out of equilibrium jsl] - ap- 
plies equally well to pure continuous media. This has been 
employed recently to transport a liquid or a solid phase in set- 
tings that do not show any macroscopic gradient The 
local asymmetry can result, for instance, from a periodic but 
asymmetric external potential. The periodic asymmetric vari- 
ation needs to be on a small length scale compared to the sys- 
tem size. 

In a first example, drops that are either placed on an asym- 
metrically structured substrate or between two such substrates 
move on average if a transverse electric field is applied peri- 
odically Jit]. A similar transport may be achieved by vibrat- 
ing the structured substrate in a tangential direction §]. A 
second example is the triggering of a large scale mean flow 
in Marangoni-Benard convection over a solid substrate with 
asymmetric grooves J^. That is an interesting example (that 
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has not yet been studied in detail theoretically) as the static 
ratchet profile of the substrate interacts with the time-periodic 
motion of the convection rolls to produce the mean flow, i.e., 
the time-periodic 'switching' is done by the system itself. The 
strength and direction of the mean flow depend on the thick- 
ness of the liquid layer and the applied temperature gradient 
across the layer. A third example are Leidenfrost drops that 
are placed on a hot surface with a ratchet-like topography. 
This induces a directed motion of the drops |0,|9t]. The effect 
is not only observed for many liquids but as well for small 
blocks of dry ice lllOll . 

In contrast to the case of particle ratchets, not many models 
exist for the ratchet-driven transport of a continuous phase. 
Reference ifTHl considers a channel flow that is induced by 
locally asymmetric periodic arrays of electrodes under an AC 
voltage ||l2il , and proposes an electro-osmotic model. 

A different concept was analysed in relation to the first ex- 
perimental example: a two-layer film confined between two 
parallel plates can be driven by periodically switching on and 
off an external potential that is of ratchet shape in space. In the 
case that the potential is an electrical one and the two plates 
form a (heterogeneous) capacitor one may call the setup a 
flashing '(electro-)wettability ratchet' iflslTTill . The concept 
of wettability may be understood in a rather general way as 
it may include any effective interactions between the liquid- 
liquid free interface and the solid walls as long as one is able 
to apply them in a time-periodic, spatially periodic (but locally 
asymmetric) manner. 

Net transport of a continuous phase may not only be caused 
by an external ratchet potential, but as well by a self -ratcheting 
effect: An enlightening example are drops that are driven up 
an incline by harmonic substrate vibrations at a finite angle 
to the substrate normal lfT5i[l6ll . A simple ratchet-like mech- 
anism has been proposed: The vibration component that is 
orthogonal to the substrate strongly modulates the hydrostatic 
pressure and thereby induces a nonlinear response in the drop 
shape. That in turn determines the strongly nonlinear drop 
mobility. As a result the drop reacts in an asymmetric manner 
to the vibration component parallel to the substrate. This sym- 
metry breaking between back and forth motion leads to the 
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observed net motion of the drop. A second experiment does 
not employ an oblique vibration but decouples the normal and 
parallel vibrations entirely iIitIi . On a horizontal substrate a 
net transport may be induced in either direction depending on 
the phase-shift and amplitude ratio of the two harmonic vibra- 
tions. Note, that in Ref. UtIi both vibrations have the same fre- 
quency. In both vibration experiments the phenomenon might 
be seen as a rocked self-ratcheting [ 1 ] of the drop as it is the 
drop itself that introduces the local time-reflection asymme- 
try in the response to the time-periodic driving of the sliding 
motion. 

The decisive element in both presented systems - the 
ratchet capacitor and the vibrated drop - is the interaction 
of the external periodic forces with capillarity and wettabil- 
ity. With other words the nonlinearity that is necessary for a 
net motion results from interface effects. One may call this 
class of ratchets, "interfacial flow driven ratchets". As they 
are dominated by interface effects they become more effec- 
tive the smaller the involved scales are. This implies that they 
are good candidates for micro- or even nano-fluidic actuators. 

In the present contribution we will restrict our attention 
to one-layer films (i.e. liquid under gas) and formulate the 
thin film model that applies to both driving types - an ex- 
ternal ratchet potential and the vibration (Section |Il|i. Then 
we present selected results for the external electrical ratchet 
potential in Section HIH Section HV] discusses drop transport 
by an oblique vibration of the substrate (cf. ifisll ). whereas 
Section |V] describes results for the recently investigated case 
where the drop moves as a reaction to a phase shift between 
decoupled horizontal and vertical vibrations Iil7i1 . Finally, we 
conclude and give an outlook in Section IvTl 



II. LONG- WAVE FILM PROFILE EVOLUTION 
EQUATION 
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The behaviour is then controlled by the interplay of the time- 
periodic external force (and its spatial modulation), wetta- 
bility and capillarity. In the limit of small surface slopes 
(small contact angles) it can be well described using an evo- 
lution equation for the film thickness profile h{x, t) that is 
derived from the momentum transport equations with ade- 
quate boundary conditions employing long-wave approxima- 
tion ifisl 191. In dimensionless form we find 

dth^ -a, {Q{h)d^ [d^^h + P{h, X, t)] + Q{h)F{t)} . 

(1) 

As mass is conserved (no evaporation is considered), the time 
derivative of the film thickness profile equals the divergence 
of a flow. The flow results from a pressure gradient and a 
lateral force and is proportional to a mobility Q{h). The pres- 
sure contains the Laplace or curvature pressure and a term 
P{h, X, t) that stands for all other pressure contributions as, 
e.g., disjoining pressure (wettability), hydrostatic pressure, 
electrostatic pressure and so on. The pressure P might de- 
pend on position and time. The lateral driving force F may 
include gravity (drop on incline), thermal or wettability gradi- 
ents and other Depending on the particular effect included it 
might as well depend on h. 




FIG. 2: (Color online) Panels (a) to (c) illustrate the working princi- 
ple of a fluidic ratchet based on a switchable wettability that causes 
dewetting-spreading cycles, (d) illustrates the spatially asymmetric 
periodic interaction profile <I>(x) responsible for the wettability pat- 
tern and (e) indicates the time dependence Q{t) of the switching in 
relation to the dewetting and spreading phases in (a) to (c). 




FIG. 1 : (color online) (a) Sketch of a film of dielectric liquid of thick- 
ness h{x) in a capacitor of gap width d and voltage U{x,t). The 
relative dielectric constants of the liquid and air are euq and Sair, re- 
spectively, (b) Sketch of a drop on a vibrating inclined substrate. The 
frequency and angle of the vibration with the substrate normal are uj 
and /3, respectively. The inclination angle of the substrate is a. 

We restrict our attention to a two-dimensional system and 
consider films/drops on lengthscales in the micrometer range. 



We are here particularly interested in three cases: 
(i) For a liquid film in an external electrical ratchet potential 
on a horizontal substrate one has 

P{h,x,t) =nvdwW + ^{t)<^{x)Ii,i{h) (2) 

and F = 0. The disjoining pressure Ilvdw comprises the ef- 
fective interactions between the liquid- gas interface and the 
substrate, i.e. the wettability properties 1 1 SL l20l I21I1 . We as- 
sume that a dielectric oil forms a film in a capacitor of gap 
width d. The oil wets the lower plate and does not wet the 
upper plate corresponding to 
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with the dimensionless Hamaker constants Ai > and Au < 
0. The dielectric film is subject to an electrical 'disjoining' 
pressure |23t] 

[Sr + (1 - er)h]^ 

that varies in space as described by $(a;) and is periodically 
switched in time as n{t). In Eq. dUi, Er denotes the ratio of 
the relative dielectric constants of the liquid, e^q, and the gas 
phase, Eair- The function $(a;) models an electric field that is 
periodic in x, but with an asymmetric profile, i.e., there exists 
no reflection symmetry in x. The employed spatial variation 
and the temporal modulation 0(<) are defined by the 
sketches in Figs|2d) and|2je), respectively. 

The introduced scales are S^rj/duli, ^^djKfa, and d for 
time t, position x, and film thickness h, respectively. Thereby, 
we have defined an electrostatic 'spreading coefficient' n^i = 
eoeiiqC^^/2rf^ where Eq is the absolute dielectric constant and 
I] is the applied voltage. 7 and r\ denote surface tension and 
dynamic viscosity of the liquid, respectively. 

For simplicity we assume Ai = —Au = v4 > 0. The di- 
mensionless Hamaker constant A is related to the dimensional 
one by A = y4.dim/67r(i'^Kei- To characterise the ratchet we in- 
troduce the flashing ratio x — '''/{T — t) and the asymmetry 
ratio (f) — \/{L — A). The flashing frequency is w = 2tt/T. 
The net transport along the substrate is measured by the mean 
flow {j). For more details see lfl3l[l4ll . 

(ii) In the second example we look at drop transport by an 
oblique vibration that determines both, P and F. The pressure 

P{h,x,t)=n{h)-Gh[l + a{t)] (5) 

contains the disjoining pressure T\-{h) = —\/h^ + l/h^ and 
the hydrostatic pressure where the time-dependence results 
from the vibration component normal to the substrate. 11 con- 
tains long-range destabilising and short-range stabilising van 
der Waals interactions i24tl . The lateral force 

F = G[a + pb{t)] (6) 

contains a constant part (force down the incline) and a time- 
modulated part (vibration component parallel to the sub- 
strate). For a harmonic oblique vibration, the substrate accel- 
eration is a{t) ~ b{t) = oq sin(ajt). Note, that the physical 
vibration angle is 0(9^13) where 9^ is the mesoscopic equilib- 
rium contact angle, i.e., a scaled angle /3 of order one corre- 
sponds to a small physical angle. 

In this case, the introduced scales are S^rj/hQii^, y^jho/n, 
and = {B/\A\)^/^ for time t, position x, and thickness h, 
respectively. A < and > are the Hamaker constants for 
the long- and the short-range part of the disjoining pressure, 
respectively. Furthermore, k — |A|/67r/ig, and G — pgho/n, 
where p denotes the density of the liquid and g the gravita- 
tional acceleration. The non-dimensional vibration period is 
T = 2tt/uj. The fixed drop volume V = L{h — hp) is deter- 
mined by the domain size L, the mean film thickness h and 
the dimensionless precursor film thickness hp = 1. The re- 
sulting transport along the substrate is measured after all tran- 
sients have decayed and the vibration-induced shape changes 



of the drop are completely periodic in time. We quantify the 
transport by the mean drop velocity {v) = Ax/T where Ax 
denotes the distance the drop moves within one period T. 

(iii) The final example is closely related to the second one. 
In the experiments of Ref. iflTil the normal and lateral substrate 
vibrations are mechanically decoupled and may have different 
amplitudes, frequencies and phases, i.e., the normal vibration 
in Eq. (|5]l is a{t) — og sin(a;at + 5), and the lateral one in 
Eg. ^ is I3b{t) = /3ao sin(ajbt). In the particular case of 
lUTIl . a horizontal substrate is used, i.e., a = and = 
LOa = u!. The parameter f3 takes the role of the ratio of the 
vibration amplitudes in the directions parallel and normal to 
the substrate. To be consistent with the long-wave approach 
taken, the physical amplitude ratio has to be small. However, 
as in (ii) (3 is the scaled ratio and therefore of 0(1). 



III. FLOW IN A RATCHET CAPACITOR 

Fig. |2] sketches an idealised electrical wettability ratchet. 
The working principle is as follows. The flat free surface of a 
film of a dielectric liquid that wets the lower wall of the capac- 
itor is stable when the electric field is switched off (Fig.|2](a)). 
When switching on the spatially inhomogeneous electric field 
at t = (see spatial profile and time dependence in Figs.|2](d) 
and|2](e), respectively) the film dewets the substrate and de- 
cays into a set of drops (we call this the 'on-phase'). This 
is due to a destabilisation of the surface by the overall elec- 
tric field and its gradients parallel to the substrate (Fig.|2](d)). 
The latter interfere with the wavelength selection in the linear 
phase of the surface instability. They do as well accelerate 
the coarsening process. The different processes during the 
on-phase can be well appreciated in Fig. |3] The qualitative 
behaviour in the on-phase of the cycle resembles dewetting of 
a liquid film on a substrate with a chemical wettability pat- 
tern llzsi I26I1 . If the on-phase is long enough (as is the case 
in Fig. O all the liquid collects in drops close to the posi- 
tions of the maximal voltage (Fig. |2](b)). Then the field is 
switched off ('off-phase') at t ~ t (Fig.|2](e)) and the drops 
spread (Fig. [2](c)), merge and become a homogeneous film 
again (Fig. [l](a)). The next cycle starts at t = T. If the pe- 
riod of the cycle is not large enough, the surface may not be- 
come entirely flat, a small modulation remains as indicated in 
Fig. He). 

When initially switching on the device there are several cy- 
cles (typically ^ 5 — 50, the number varies with the employed 
parameters) that show transient behaviour However, the ini- 
tial transients die out and the evolution of the film profile dur- 
ing one cycle is exactly time-periodic. Note, that there exist 
small parameter ranges where this is not the case due to reso- 
nance phenomena and/or multistability of various film states. 
However, this shall not concern us here. 

Although the evolution of the film profile is exactly time- 
periodic, the process re-distributes the liquid within the film. 
This results in a net transport that depends in a non-trivial way 
on the various control parameters. As an example we give in 
Fig- SI a log-log plot of the dependence of the mean flux on 
the flashing frequency. For the used flashing ratio and ratchet 
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FIG. 3: Shown is a space-time plot of the typical evolution of the 
film thickness profile of a film of dielectric liquid in a heterogeneous 
capacitor (one spatial period shown) during one flashing period. For 
capacitor geometry and flashing cycle see Fig. |2l Parameters are 
h = 0.5, * = 0.5, L = 32, fli = 5, T = 5000, x = 1, ^ = 
0.001, Er = 2.5. The starting time is well after initial transients have 
decayed. 




FIG. 4: (Color online) Shown is the net transport as measured by 
the mean flux (j) in its dependence on the flashing frequency co for 
various film thicknesses as indicated in the legend. The remaining 
parameters are as in Fig. [3] The dotted straight line corresponds to a 
power law of exponent one. 



geometry the flux is always positive and approaches zero in 
the zero frequency Umit as wefl as in the high frequency Hmit. 
For large frequencies the fluid does not have enough time to 
dewet in the on-phase or to spread in the off-phase. At small 
frequencies both processes reach the respective equilibrium 
well before the next switching, i.e. most time is spend waiting. 
As a result the mean flow increases proportionally with fre- 
quency when increasing the frequency from zero (cf. Fig.|4|. 
The flux reaches a maximum at w = lO^'^ . . . 10~^ before it 
decreases again. The decrease seems to follow a power law 
as well. It is, however, not universal: the exponent increases 
with increasing mean film thickness. The decaying part shows 
non-monotonic step-like behaviour that is more pronounced 
for smaller film thicknesses. This results from that fact that 
films of smaller thicknesses dewet in a homogeneous electric 
field with a wavelength well below the spatial period of the 
ratchet. The more wavelengths fit into a ratchet period the 
more coarsening steps can take place in the on-phase if the 



frequency is small enough. Each of the 'steps' in the curves 
of Fig.lHis related to a coarsening step that does not take place 
above the frequency at the step. At smaller mean film thick- 
nesses the maximum resembles a plateau with a flux that does 
nearly not change between ui = 10^^ and 10^^. The maxi- 
mal flux increases with increasing mean film thickness. The 
particularly interesting non-monotonic behaviour close to the 
flow maximum (here well visible for h = 0.5) results from 
the interaction of the dynamics of the 'last' coarsening step 
and the switching cycle. 

Second, Fig. |5] shows the dependence of the mean flux on 
the flashing ratio. When it is low, i.e., the ratchet is most of the 
time in the off -phase, the flux is very low. This is because the 
liquid remains nearly homogeneously spread out as it has no 
time to assemble at the points of maximal voltage. By increas- 
ing the time of the on-phase the flux increases in a non-trivial 
manner until it reaches a maximum at x ~ 1- Depending on 
the film thickness one may find almost a plateau, i.e. a range 
of flashing ratios with nearly constant flux. For high flashing 
ratios the flux decreases again, since the ratchet is most of the 
time in the on-phase. That implies electrically formed drops 
have not enough time to spread out. The decrease follows a 
power law with an exponent of about —4/3. 
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FIG. 5: (Color online) Shown is the net transport as measured by the 
mean flux (j) in its dependence on the flashing ratio x for various 
film thicknesses as indicated in the legend. T = 500 and the remain- 
ing parameters are as in Fig. [3] The dotted straight line corresponds 
to a power law of exponent -1.3. 



Further calculations (not shown) indicate that the mean flux 
increases monotonically with increasing spatial asymmetry 
ratio (f> and the amplitude 5* of the ratchet potential. Note 
that as expected the net transport is zero at = 1, i.e., for a 
symmetric potential. For further results for the presented one- 
layer geometry see Refs. ifisi [fl . Exchanging the air layer 
by a second dielectric liquid one is able to transport the two 
liquids into opposite direction. There, depending on the ratios 
of the viscosities and relative dielectric constants one may as 
well find a flux reversal. This is further discussed in Ref. lfl4ll . 
With this we close our present discussion of case (i), i.e., the 
study of the transport of a dielectric liquid by a ratchet ca- 
pacitor Our results are typical for many externally imposed 
ratchets. The two cases discussed next are both related to net 
transport through substrate vibrations. Note that they can be 



5 



seen as 'intrinsic' or self-inflicted ratchets. 



IV. DROP TRANSPORT BY OBLIQUE SUBSTRATE 
VIBRATION 




FIG. 6: (color online) Space-time plot illustrating the evolution of 
the profile of a drop on a obliquely vibrated substrate during one 
vibration period. Each cycle results in the net motion of the drop to 
the left. The starting time is well after initial transients have decayed. 
Note that only part of the domain L is shown. Parameters are V = 
192, G = 0.001, /3 = 0.1, Q = 0, T = 400 (a; ^ 0.017), L = 128, 
ao = 10. 

As the first case of vibration-induced transport we discuss 
case (ii) introduced in Section HI]- a drop transported by an 
oblique substrate vibration as sketched in Fig. [Hb) and ex- 
perimentally observed in Ref. [15]. The typical behaviour of 
a drop during one vibration cycle as obtained in a time sim- 
ulation of Eq. ([T) is given in Fig.|6] One observes that the 
drop undergoes changes of shape and moves back and forth. 
Such simulations in time may be employed to analyse the net 
transport over a wide frequency range. In addition one may 
use continuation techniques ll27l - i29ll to study the behaviour in 
the low frequency limit, i.e., for a slowly vibrating substrate. 
In this limit the intrinsic timescale of the drop dynamics is 
much smaller than the vibration period T and the drop moves 
in a quasi-stationary manner This means that drop shape and 
velocity at each instant during the vibration cycle correspond 
to the ones of a stationary moving drop at the corresponding 
constant force. They are parametrized by a{t). Averaging sta- 
tionary drop velocities v{a{t)) over one vibration period gives 
the low frequency limit of {v). 

Fig.|7]presents results for the low frequency limit. Panel (a) 
shows profiles of stationary moving drops at various phases 
of the cycle on a horizontal substrate. Panels (b) and (c) give 
drop velocity and the dynamic mesoscopic contact angles, re- 
spectively, over one vibration cycle. The angles are shown 
for the left (6*;) and right {9^) contact line. In the given setup, 
temporal modulations of the right contact angle are larger than 
those of the left one. The deviation of the velocity of the slid- 
ing drop from a harmonic modulation [dotted line in Fig.|3b)] 
is a measure of the net motion of the drop. 

The overall behaviour of the drop over one vibration cycle 
is very similar in Fig. |6] (finite, but small frequency) and Fig.|7] 
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FIG. 7: (color online) Shown are for an obliquely vibrated horizon- 
tal substrate in the low frequency limit, (a) several drop shapes at 
times as given in the legends, (b) the momentaneous velocity, and 
(c) the left {6i) and right {9r) mesoscopic dynamic contact angle dur- 
ing one vibration cycle. For the used /? = 0.3 the drop moves with 
[v) ~ —0.01, i.e., to the left. The blue dotted line in (b) indicates a 
harmonic variation of zero net flow. Remaining parameters are as in 

Fig.m 



(low frequency limit): during the first half of the cycle (t < 
T/2) the drop is flattened and moves to the right. Note that 
this happens when the substrate is accelerated upwards and to 
the left. In the second half of the cycle the drop becomes taller 
and less wide while it slides to the left {t > T/2, substrate ac- 
celeration is downward and to the right). After one period the 
drop has moved a small net distance to the left. For the param- 
eters of Fig.|6]it takes about 100 vibration cycles to move the 
drop by its own length. As detailed below in the conclusion, 
this does well correspond to the available experimental results 

iH. 

Based on the described findings one is able to understand 
the mechanism that leads to the net motion of the drop: The 
component of the oblique vibration that is normal to the sub- 
strate strongly modulates the hydrostatic pressure and there- 
fore provokes a nonlinear response in the drop shape, i.e., 
when the substrate accellerates upwards [downwards] it com- 
presses [decompresses] the drop. That in turn determines the 
strongly nonlinear drop mobility Q{h) in Eq. ([U and is there- 
fore responsible for an asymmetric response to the back and 
forth forcing that results from the parallel vibration compo- 
nent. As a result of this nonlinear coupling of the effects of 
normal and parallel vibration component one obtains an an- 
harmonic response of the drop to the harmonic but oblique 
vibration. 

The examples given in Figs. |6] and [T] where for a rather 
small vibration angle /3. Therefore the change in drop shape 
was mainly due to the modulation of the hydrostatic pressure. 
For larger /3 the drop additionally changes its shape in re- 
sponse to the parallel vibration component. An examples is 
given in Fig. [8] The characteristic backwards shoulders are 
known from strongly driven drops on homogeneous substrate 
ll29l[30ll . For even larger /3 the process may become very com- 
plicated as the drop can during a single cycle undergo several 
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FIG. 8: (color online) Shown are several drop shapes during one 
vibration cycle for an obliquely vibrated horizontal substrate alV — 
192, G = 0.001, ao = 10, a = 0, r = 400, L = 128 and a 
relatively large vibration angle 13 = 1. 
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FIG. 9: Scaled mean velocity depending on the vibration frequency 
uj for drops on a horizontal substrate. Panel (a) gives the master 
curve for the scaled mean velocity {v) //Soq for sets (ao, /3) as given 
in the legend. The horizontal dashed line indicates the result in the 
low frequency limit. Remaining parameters are as in Fig. [6] Panel 
(b) gives the scaled mean velocity {v)/V^'^''^ for drops of different 



volume as given in the legend. Remaining parameters are ao 
/3 = 0.1, Q = 0.0, G = 0.001. 
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FIG. 10: Drop profiles at various times of the cycle as indicated in 
upper right corner of each plot for various vibration periods T as 
indicated in the legend. T = oo indicates the low frequency limit. 
Remaining parameters are V = 384, L = 256, a = 0, /? = 0.1, 
ao = 10, G = 0.001. 
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FIG. 1 1 : Variation of the maximal drop height during one vibration 
period for different vibration periods as indicated in the legend. Re- 
maining parameters are as in Fig.llOl 



morphological changes. It can transform between a spherical 
cap-like drop via a drop with a backwards shoulder to a fi- 
nite film with a capillary rim [30]. We do here, however, not 
investigate this regime further. 

For smaller vibration angles (3 morphological transitions do 
not occur (at reasonable accellerations ao)- This allows for 
'universal' behaviour, i.e., one finds that scaling laws that are 
discovered in the low frequency regime hold in part in the en- 
tire frequency range studied. In particular, one finds in the low 
frequency regime a scaling {v) ~ fSa^L^-^'^. Fig.|9ja) shows 
that for drops of identical volume the scaling (v) ^ I3a^ does 
hold very well as it is possible to 'collaps' curves for a range 
of parameter pairs (oq, /3) on a single master curve. The scal- 
ing with volume does, however, not hold [see Fig. |9jb)]. In- 
terestingly, one finds a flux reversal at high frequencies, i.e., 
above a critical frequency ujc the drops move on average to the 
right. Fig. |9jb) further indicates that for larger drop volume 
the reversal is more pronounced and occurs at lower frequen- 
cies LOc- We will end this section with an investigation of the 



To this end we show in Fig.[TO]a superposition of drop pro- 
files obtained at frequencies above and below loc at selected 
times during the vibration cycle. The profiles in the low fre- 
quency limit are per definition ideally in-phase with the sub- 
strate vibration. At T = 400 (a; = 0.016 < « 0.05, 
cf. Fig.|9|b)) the profiles near the times of maximal acceller- 
ation (i.e., t sa T/A and t « 3T/4) strongly resemble the 
ones in the low frequency limit, although the variation in their 
shape is already smaller For T = 1000 they are nearly in- 
distinguishable (not shown). However, around t = Q and 
T = T/2 where the accelleration changes sign the conver- 
gence to the low frequency limit is much slower as there even 
small phase lags are very important (cf. Fig.fTTI). The overall 
behaviour can be well appreciated in Fig.[TT]where the depen- 
dence of the maximal drop height on time is shown for one 
vibration cycle for several different periods. With decreas- 
ing period we observe a continuous increase of the phase shift 
w.rt. the curve for the low frequency limit. The curves of 
finite frequencies lag behind the one for zero frequency and 
show a less pronounced variation of the drop heights. 
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FIG. 12: (color online) Shown is the mean velocity as a function of 
the phase shift 5 for various amplitude ratios j3 as given in the legend. 
Thereby, we use a relatively large period of T = 200 and give the 
scaled velocity For j3 < 0.2 all curves 'collapse' onto a 

master curve. The remaining parameters are V = 192, G — 0.001, 
ao = 10, a = 0, L = 128, ao = 10. 



At T = 70 (w = 0.090 > ojc) the drop profiles in Fig.[TO] 
vary much less over time than for the lower frequencies and 
Fig, nn shows a rather large phase shift (larger than 7r/3) as 
compared to the other curves. When increasing the frequency, 
due to the phase shift and nonlinear mobility eventually the 
average response to the backwards force becomes smaller than 
the one to the forward force and the drop shows a reversal of 
the net motion. The net motion to the right becomes maxi- 
mal at an optimal frequency. The influence of the phase lag 
still exists at higher frequencies, however, as the changes in 
the profile shape become smaller with increasing frequency, 
the mean velocity decreases again and approaches zero in the 
limit of large frequencies. This is as well the reason why the 
absolute values of the velocities that can be reached in the 
flow reversal regime are much smaller than the ones in the 
low frequency limit. Note, finally, that in the case of larger 
vibration angles the picture becomes more complicated as the 
phase lag interacts with the morphological changes. As a re- 
sult the dependency of mean velocity on frequency might have 
more than one maximum. 



V. DROP TRANSPORT BY DECOUPLED HORIZONTAL 
AND VERTICAL SUBSTRATE VIBRATIONS 

As a second example of vibration-induced transport we dis- 
cuss case (iii) introduced in Section|II]- a drop transported on 
a horizontal substrate by decoupled normal and parallel sub- 
strate vibrations as experimentally observed in Ref. UtIi . We 
focus on normal and parallel vibrations of identical frequency 
that have a phase shift and possibly different amplitudes. Fol- 
lowing the discussion of the underlyig mechanism in case (ii) 
it is no surprise that a net motion is found with the present 
simple long-wave model as well for case (iii). 

Fig. [T2I shows the dependence of the scaled mean velocity 
on the phase shift for different amplitude ratios /3 at a moder- 
ately large vibration period of T = 200 (lu = 0.031). Very 
similar behaviour is found for other periods. Interestingly, the 
dependencies of the scaled velocity {v)//3 for different /3 fall 



FIG. 13: (color online) Shown is the mean velocity as a function 
of the phase shift 5 for various periods as given in the legend. The 
remaining parameters are V — 192, G — 0.001, ao = 10, P — 0.1, 
a = 0, L = 128, ao = 10. 



on a single master curve for /3 < 0.2. This implies that for 
/? < 0.2 particularly important values of the phase shift do not 
depend on the amplitude ratio. In particular, we find that flow 
reversal occurs at S = 3/47r and 6 = 7 /An. The drop moves 
fastest to the left [right] at S = tt/A [S = 5/4tt]. The master 
curve shows the symmetry {6 6 + it,{v)/13 — 
Note that the symmetry holds for all considered amplitude ra- 
tios and periods (cf. Fig. \T% . as can be expected from the 
set-up of the problem. 

However, the scaling found for small /3 does not hold for 
larger ^ (see Fig. [121). For 0.2 < /? < 1.0 the phase shift 
that results in maximal net motion and the maximal net ve- 
locity {v) both increase with increasing /3. For /3 > 1.0 both 
decrease again (not shown). This implies that the amplitude 
ratio that maximises net transport is /3 1. There the phase 
shift resulting in maximal net transport to the left is slightly 
smaller than S = Tr/2. 

The change of the dependency of transport on phase shift 
with changing period is shown in Fig. [13] In accordance with 
results in case (ii) one finds that the velocities become mono- 
tonically smaller (for all phase shifts) with increasing frequen- 
cies (decreasing periods) as the drops are less able to follow 
the substrate vibrations. The phase shift where the maximal 
mean velocity occurs becomes smaller with decreasing fre- 
quencies. In the limit of low frequency the fastest transport to 
the left and right are found atS — and S — tt, respectively. 

A comparison with the experiments in Ref. ITtIi (in partic- 
ular their Fig. 3) shows that there the phase shift value where 
the maximal mean velocity is found and the value of the max- 
imal mean velocity both increase with increasing ratio of the 
amplitudes of parallel and normal vibration. However, as fur- 
ther explained in the conclusion, a direct quantitative compar- 
ison is not possible. The experimental results show roughly 
the symmetry {6 ^ S + tt, (v) — ?> —(«))• Small deviations 
are explained by the presence of defects on the substrate UtIi 
(see, in particular, the inset of their Fig. 2). A striking differ- 
ence between the results presented here and Fig. 3(a) of UtIi 
is that in the latter the mean velocity depends in a strongly 
non-harmonic way on the phase shift. In the simulations we 
have not seen the experimentally found double peak structure 
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around the maximal velocity. This point needs further investi- 
gation. 



VI. CONCLUSIONS 

We have explored the usage of ratchet mechanisms to trans- 
port a continuous phase in micro-fluidic settings involving a 
free liquid-gas interface and contact lines, i.e., under the in- 
fluence of capillarity and wettability. In particular, we have 
studied on the one hand the transport of a dielectric liquid in 
a capacitor with an asymmetrically spatially modulated elec- 
trical field that is periodically switched on and off. In this 
case the ratchet-like potential is imposed externally. On the 
other hand we have investigated drops on a solid substrate that 
show net motion under the influence of different types of har- 
monic substrate vibrations. Analysing the underlying mecha- 
nism we have found that the component of the harmonic vibra- 
tion that is orthogonal to the substrate induces a nonlinear (an- 
harmonic) response in the drop shape. The latter determines 
the strongly nonlinear drop mobility what results in an asym- 
metric response of the drop to the vibration component that 
is parallel to the substrate. The induced symmetry breaking 
between forward and backward motion during the different 
phases of the vibration results in the observed net motion of 
the drop. We have argued that the phenomenon might be seen 
as a rocked self-ratcheting as the drop itself introduces the 
local time-reflection asymmetry in the response to the time- 
periodic driving of the sliding motion. We employ the term 
"rocked" as the net motion is strongest in the low frequency 
limit as is typical for rocked particle ratchets [1]. Note as well 
that in the context of a drop on an incline the component of 
the vibration parallel to the substrate can be seen as a periodic 
rocking of the substrate. 

For the vibrated drops our focus has been on two experi- 
mentally realised cases ITsUitIi : an oblique substrate vibra- 
tion and decoupled normal and parallel substrate vibrations 
of identical frequency but different amplitudes and/or phase 
shifts. As in the latter case the use of different amplitudes 
without phase shift corresponds to the case of the oblique vi- 
bration we have mainly analysed the dependence of the drop 
motion on phase shift. 

The decisive element in both presented systems - the 
ratchet capacitor and the self-ratcheting vibrated drop - is 
the interaction of the external periodic forces with capillar- 
ity and wettability, i.e., the nonlinearity that is necessary for 
a net motion results from interface effects. One may call this 
class of ratchets, "interfacial flow driven ratchets". As they 
are dominated by interface effects they are very effective on 
small scales. This implies that they are good candidates for 
micro- or even nano-fluidic actuators. 

One needs to keep in mind that all models analysed in the 
present contribution are obtained employing a long-wave or 
lubrication approximation ifislfigll . That means they are rela- 
tively simple, but fully dynamical, highly nonlinear models 
based on a minimal set of ingredients, i.e.. Stokes flow in 
the lubrication approximation, (electro-)capillarity, and wet- 
tability. They allow to study the net transport in the ratchet 



capacitor and the net motion of the vibrated droplet over a 
wide range of parameters, to understand the underlying mech- 
anisms and to discuss transitions in the qualitative system be- 
haviour. 

However, when comparing such lubrication models with 
particular experiments one always has to keep in mind that 
the formal range of applicability of the lubrication approxima- 
tion is limited and does often not overlap with the parameter 
ranges where experiments work best. The latter are often per- 
formed with liquid-substrate combinations that lead to equi- 
librium (and dynamic) contact angles that are not small. In 
our case (ii) one uses as well large substrate inclinations and 
large angles between vibration direction and substrate normal. 
For a direct quantitative comparison all those angles have to 
be small. Nevertheless, lubrication models are extremely suc- 
cessful in explaining intriguing effects observed in a wide 
range of experiments involving capillarity and wettability. A 
good example for this are studies of morphological changes 
in sliding drops where some experiments are done for liquids 
with static contact angles above 7r/4 using substrate inclina- 
tions varying from zero to tt/ 2 Olll . Lubrication theory well 
explains the effects (see, e.g., ||32I1 and the discussion in |[33il ). 

For the vibrated drop in Ref. ifTsIl a liquid is used that has 
an equilibrium contact angle of about 7r/3, the substrate incli- 
nation is 7r/4, and the vibration is applied at an angle of 7r/4 
with the substrate normal, i.e., the experiment is well outside 
the formal range of applicability of the lubrication approxima- 
tion. Therefore, a direct quantitative comparison of mean ve- 
locities could only be done through some 'up-scaling' proce- 
dure that lifts drop volume, contact angles and vibration angle 
from the lubrication theory values to the experimental ones. 
This might be done (and is sometimes done). However, in our 
opinion the process has arbitrary elements, and we prefer not 
to employ it. 

The theory may, however, be related to the experiment in a 
semi-quantitative way by comparing, for instance, mean ve- 
locity in terms of drop size or typical timescales. In Ref. lITsIl 
the employed frequencies / range from 25 Hz to 120 Hz. At 
an acceleration of 15 g and / = 60 Hz one finds mean veloci- 
ties of about Imm/s, i.e., the drop moves less than 20 fim per 
cycle. For a drop of 3 mm length this implies that it needs 
about 150 cycles to move by its own length. The situation is 
very similar in Ref. iflTll where for the parameters depicted in 
their Figs. 1 and 2 it takes about 50 vibration cycles to move 
the drop by its own length. Note, that their Fig. 1 does not 
show snapshots taken during a single cycle, but selected pic- 
tures from 70 cycles. 

In terms of the number of cycles needed to move the drop 
by its own length our results are very close to the experimental 
findings: At a period of T = 100 we have a scaled velocity of 
about 0.1 X 10~3 (our Fig.^a)); for oq = 10 and (3 = 0.2 
this gives (w) = 2 x 10^^ and an advancement per cycle of 
Ax = 0.2. For a drop like the one in our Fig. |6] this implies 
that it takes about 190 cycles to move it by its own length. 
Employing instead the accelerations of 15.g as in ifTsIl it takes 
about 90 cycles. This is a fair agreement between experiment 
and lubrication theory. One may as well discuss the time scale 
to employed in the non-dimensionalisation and use it to com- 
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pare the relevant vibration frequencies in the model and in the 
experiment. 

Finally, we would like to discuss the relation of the dis- 
cussed ratchet mechanisms to transport a continuous phase 
to 'particle ratchets'. We do this by pointing out some for- 
mal analogies between the employed thin film equation and a 
Fokker-Planck equation for transport of discrete objects in a 
particle ratchet. 

The film thickness evolution equation 

dth = ~d, {Qih)d^ [d.,.h + P{h, .T, t)] + Q[h)F{t)} . 

(7) 

shows similarities to a Fokker-Planck equation for interact- 
ing particles in an external ratchet potential ll34l[35ll . Its non- 
dimensional form is 

dtW = -d, {Wd,P{W, X, t) + WFit)} (8) 

with 

P{W,x,t) = -T{t)\ogW -gW ~U{x,t) (9) 

In the particular case, it is an evolution equation for the 
one-particle distribution function W{x,t), where U{x,t) is 
a flashing ratchet potential (note that the one used in 1341 13511 
is static) 'rocked' periodically by F{t). The function T{t) 
stands for a periodic modulation of the non-dimensional tem- 
perature, and g is a non-dimensional interaction parameter. 
One notes that both, the first order drift or transport term and 
the second order diffusion term are modulated in time. Note 



that U{x, t) could as well be incorporated into the drift term 
instead of the diffusion term. 

However, differing from ||34[ Issll our cases (ii) and (iii) 
do not incorporate an external potential. The absence of a 
spatial ratchet potential is compensated by strongly nonlin- 
ear prefactors of the diffusion and the drift term that result in 
mean flow even for harmonic (i.e., time-symmetric driving). 
In the context of 'particle ratchets' related systems are studied 
ll36[ [37I1 where a mean flux is created without spatial ratchet 
potential but through the response of two different (weakly) 
coupled degrees of freedom to time-periodic driving. The de- 
grees of freedom are related to two different interacting par- 
ticle species in ll36ll . and to pancake vortices and Josephson 
vortices in layered superconductors in iIstIi . However, in con- 
trast to our case in ||36[ [stIi the driving is time-asymmetric. 
Time-symmetric driving can as well induce mean transport in 
a 'particle ratchet' without spatially varying potential if the 
coupling of the different degrees of freedom is sufficiently 
nonlinear. The only example we are aware of is the "vortex 
diode" discussed in Issll . Note finally, that similar 4th order 
terms as the capillarity term in Eq. ^ may be included in 
higher order Fokker-Planck equations. 
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